Bi-defects of Nematic Surfactant Bilayers 
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We consider the effects of the coupling between the orientational order of the two monolayers in 
flat nematic bilayers. We show that the presence of a topological defect on one bilayer generates a 
nontrivial orientational texture on both monolayers. Therefore, one cannot consider isolated defects 
on one monolayer, but rather associated pairs of defects on either monolayer, which we call bi-defects. 
Bi-defects generally produce walls, such that the textures of the two monolayers are identical outside 
the walls, and different in their interior. We suggest some experimental conditions in which these 
structures could be observed. 
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Nematic liquid crystals are fluid phases possessing a 
long-range orientational order j^]. Ordinary nematics, 
in the three-dimensional (3D) space, consist of rodlike 
molecules orienting parallel to some unit vector n, called 
the "director" . Since nematics bear no polar order, n 
and — n represent the same orientational state. Nematics 
exhibit striking (line or point) topological defects [|l],D. 
The orientational order is continuous outside the defect, 
but exhibits on it a singularity which cannot be removed 
by continuous deformations. 

Although several almost 2D nematic systems have been 
investigated, like thin nematic cells [|j and wetting lay- 
ers m, there are few examples of real 2D nematics, e.g., 
rods suspended on the surface of aqueous solutions 
(Actually, 2D systems can only exhibit guasi-long-range 
order, but this distinction is blurred for usual system 
sizes.) Very recently, it has been shown that amphiphilic 
bilayers made of dimeric surfactants (gemini) sponta- 
neously form very long tubules of mesoscopic radius 
this conformation can be theoretically explained by in- 
troducing a coupling between the surface curvature and 
two independent monolayer nematic orders j7j. A num- 
ber of independent arguments support the existence of 
nematic order in these membranes Q. 

In this Letter, we investigate the behavior of discli- 
nation defects in such nematic bilayers. For simplicity, 
we restrict our attention to planar bilayers, which could 
be produced by osmotically blowing up the tubes, or by 
patch-clamping techniques. We find radically new fea- 
tures due to the coupling of the nematic order between 
the two monolayers. Even if a disinclination is present on 
only one layer, the coupling generates a nontrivial tex- 
ture on the opposite one: this texture must be considered 
as a "defect" even in the absence of a singularity. We are 
thus led to consider pairs of associated defects on the bi- 
layers, one of which can be virtual (of zero strength). We 
call these structures hi- dejects. 

We show that the two interacting nematic monolayers 
can be mapped on two independent, "virtual", 2D ne- 
matic monolayers, one subject to an external orienting 



field and the other free. In the former, defects generate 
orientational walls i.e., ribbons where the director 
turns by TT on a finite length. Consequently, bi-defects 
generally produce walls that reach the boundary of the 
sample: the textures of the two monolayers are identi- 
cal outside the walls and different in their interior. The 
bi-defect energy is dominated by the walls, and scales 
therefore linearly with the sample size (rather than log- 
arithmically) . 

We denote by m and n the directors of the upper and 
lower monolayer, respectively. Within the one Frank con- 
stant approximation, the nematic free energy of the bi- 
layer can be written as 

F = i Jd^r^^K |Vm|^ + K |Vn|^ - A (m • n)^} , (1) 

where, e.g., |Vnp = diUj diUj and summation on re- 
peated indices is understood. To be definite, we suppose 
A > 0. This is no restriction, since there is always the 
freedom to redefine n by a 7r/2 rotation, which effectively 
changes the sign if the interaction term in Eq. (^ . Let us 
call 6*+ (resp. 9-) the polar angle of m (resp. n) relative to 
an arbitrary direction. Setting 9± — ^(^±■0), we obtain 
(up to an irrelevant additive constant) F — ^(-Fb + -Fa), 
with 

Fa = j d'r I {V^f ■ (2a) 

Fx = Jd^r^Y ^^'^^^ + ^ ''''''^} • ^^^^ 

Equation ( pa| ) describes a free nematic, while Eq. ( pb| ) 
describes a nematic subject to a uniform field directed 
along the ■0 = axis The Euler-La grange equation 
deriving from ( p^ ) is a sine-Gordon equation: 

C2v2(20) = sin(20), (3) 

where — K/{2\). The length ^ is the analog of the 
magnetic coherence length of ordinary nematics Q . The 
corresponding equation for (f> is simply = 0. 
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A topological defect of strength p, located at the ori- 
gin, is described in polar coordinates by solutions of the 
Euler-Lagrange equations of the form 

^ir,9)^p0 + Mr,O), (4) 

(or the analog for ip), where p is a half-integer, and 
(j)c{r,6) is a continuous function. Indeed, the director 
turns by 2p7r in any circuit around the origin. In the 
nematic under field, minimization of the energy requires 
that = kn (where k is an integer) over most of the 
sample. Therefore, all nonuniformity is confined within 
"soliton" walls of thickness ~ 5^, crossing which ip ro- 
tates by ±7r P,|l|. Thus, a defect of strength p radiates a 
"star" of 2\p\ walls. Within a region of size ~ ^ around 
the defect the texture is similar to that without field. 

In mean field, the energy of a defect of strength p in a 
free nematic is equal to TrKp\n{L/a) [||, where L is the 
linear size of the sample and a the radius of a core inside 
which the nematic order is destroyed. The interaction 
energy of two defects of strength pi and p2 is given by 
—2TrKpip2ln{d/a), where d is the distance between the 
defects Q . For the nematic under field, the defect energy 
is dominated by the energy of the walls, which is equal 
to 2K/^ per unit length 0. 

Since Eq. (|^) is linear in the defect strength, a bi-defect 
\p,q], i.e., the superposition of a defect with a strength 
p in the upper monolayer and a strength q in the lower 
one, is equivalent to a pair of defects of strength p + q in 
the free nematic (described by 0) and of strength p — q 
in the nematic under field (described by ip): 
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We call p-\-q the free strength and p — q the field strength 
of the bi-defect. It follows from our decomposition that 
a bi-defect of free strength £ and field strength m obeys 
the relations 

^±({m}) =^(^o(^)±0a(™)); (6a) 

^({m}) = \{m)+F,{m)). (6b) 

In this equation, 9q{1) is the texture of a defect of 
strength I 'ma free nematic, 9\{m) the texture of a defect 
of strength m in a nematic under field, and -Fb(^), F\{in) 
the corresponding energies. In particular, 9^{[p,Q\) = 
9^ {{p,p}) — i {9o{p) — 9\{p)), and therefore there is a 
nontrivial texture even in the lower monolayer of a [p, 0] 
bi-defect, where there is no singularity. 

By applying these rules, one can build up the textures 
corresponding to different bi-defects. Figure |^(a) shows 
the texture of a [^,0] bi-defect. The full (resp. dashed) 
lines are the field lines of the upper (resp. lower) 
monolayer, and the wall boundary is indicated by the 




(a) (b) 

FIG. 1. (a) Field lines of a [i, 0] bi-defect. The textures of 
the two monolayers coincide outside the wall, (b) Level lines 
of the corresponding {5,5} bi-defect. 




(a) (b) 
FIG. 2. (a) Field lines of a [1,0] bi-defect. (b) Level lines 
of the corresponding {1, 1} bi-defect. 

dotted line. The corresponding ^} texture is shown 
in Fig. 0(b): the bold lines are the level lines for the free 
nematic, and the thin ones for the field nematic. Figure^ 
shows the analog texture for a [1, 0] bi-defect. 

In crossing a wall, both 9+ and 9^ turn by ±7r/2. The 
actual thickness of the walls is ~ 5^, as one finds by 
integrating Eq. (^). Besides, the walls probably have a 
persistence length which is several times their thick- 
ness. Since they are interfaces in two dimensions, they 
fluctuate widely: their lateral excursion Au over a length 
L is given by 




where T is the temperature, measured in energy units. 
They perform therefore a random walk, but their angular 
fluctuation Aa ~ (T/4i4:)i/2(f /$p)i/2 is smah, since we 
expect K to be of order a few T in a nematically ordered 
phase. The fluctuations of the wall decrease the effective 
line tension by a negligible amount. 

The walls issuing from bi-defects can recombine. Since 
a defect of strength I under field generates 2\i\ walls, a 
[p, q] bi-defect generates 

n = 2\p-q\ (8) 
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(c) 

FIG. 3. (a) Field lines of a [1, 0] & [0, 1] pair of bi-defects. 
(b) Corresponding scheme of the wall connections, (c) Possi- 
ble wall connections between three bi-defects. 



walls. Now, if there are two bi-defects, of strengths [p, q] 
& [p' , q'] respectively, the total field strength equals p — 
q + p' — q' and the number of walls that reach infinity 
is then 2\p ^ q + p' — q'\. If this number is smaller than 
2|p ~ ll + 2|p' — q'\, some walls must recombine. This 
happens if {p — q){p' — q') < 0. Therefore, we can assign 
an arrow to each wall, pointing outward from the bi- 
defect if (p — > and toward it otherwise: walls with 
matching arrows can recombine. We show in Fig. |^(a) 
the field lines of a bi-defect pair [1, 0] & [0, 1]. The two 
pairs of walls combine, connecting the two bi-defects, as 
shown in Fig. ||(b). 

The interaction energy of a [1, 0] & [0, 1] bi-defect pair, 
equivalent toa{l&l,l&— 1} system, can be estimated 
using Eq. (6b). The first contribution, ^-Fbj is one half of 
the energy of a pair of defects of strength 1 in a free ne- 
matic, i.e., — 7r/Cln(d/a), where d is the distance between 
the defects. The second contribution, \Fx, is one half of 



the energy of the texture under field of a pair of defects 
of strength 1 and —1. When d ^ ^ it is dominated by 
the two walls which connect the defects, and is therefore 
~ AK{d/£). When d ^ f we can distinguish a region of 
size « ^ where the texture is similar to that without field, 
and an exterior region where "0 is exponentially close to 
kiT. The corresponding energy ([2b|) contains two contri- 
butions: the elastic energy tt/C ln((i/a) and the potential 
energy, which is estimated by integrating ^Xsvcl' il) for 
the free texture on a disk of radius « ^. One obtains 
lK{d/£,f[\n{S,/d) -f i]. Summing up ^Fq and \Fx we 
obtain 




for d < f ; 
for d > ^. 



(9) 



The two bi-defects are therefore attracted by a force 
which is almost constant at large separation, and van- 
ishes roughly linearly with d when d <C . Indeed, when 



the two bi-defects sit on top of each other, they form a 
[1, 1] bi-defect which optimizes both the coupling and the 
elastic energies. 

A [1, 0] & [0, —1] bi-defect pair, equivalent to a {1 & —1, 
1 & 1} system, generates two walls which wander to the 
boundary of the sample. The elastic energy, calculated 
as previously, is given by 



int 



ttK In - + Fwaiis, for d > ^. 
a 



(10) 



The energy in the first line is simply the integral of the 
iAsin^i/' term. (The free and field elastic energies com- 
pensate as previously.) There is also a contribution due 
to the walls, but it does not depend on d. The first term 
in the second line represents the logarithmic attraction 
of the defects in the free nematic. -Fwaiis is the contri- 
bution from the walls of the nematic under field. It will 
depend in general on the way the walls reach the sample 
boundary. Let us consider, e.g., the case in which the 
sample is a ribbon of width 2L, with the two bi-defects 
in the middle, each sending a wall to the opposite sides 
of the ribbon. Each wall of length L wanders within a 
rectangular region of width Am given by Eq. (|). Thus, 
if d > Au, -Fwaiis is independent of d, whereas, if d < Au, 
there is a Helfrich-like repulsion between the walls: 



-^walls 



K d2 ■ 



(11) 



Therefore the interaction is repulsive for d <tc ^, and 
is otherwise a combination of repulsive and attractive 
forces, which identify an equilibrium distance 



J-cq 



(12) 



Let us now consider a collection of bi-defects [pi,qi] 
placed in a region of size R inside a sample of size L ^ R. 
Since the total field strength is given by X^iPi ~ J2i^i^ 
there are 



(13) 



walls going to the boundary. Since the total number of 
walls issuing from the defects is \pi — qi\, there are 



m=y: 



(14) 



walls linking two bi-defects, that remain confined within 
R. Therefore the dominant energy, which arises from the 
walls, scales as 



F ^ N + M Kj. 



(15) 
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In order to minimize its energy, the system will first at- 
tempt to bring N to zero, e.g., by nucleating defects on 
the boundary, so as to equalize the total strengths of the 
defects in the upper and lower monolayer. The follow- 
ing step will be to bring together bi-defects having field 
strengths of opposite sign, in order to reduce to « ^ the 
total wall length. The bi-defects can then recombine. 

In unconstrained membranes, there are numerous and 
subtle effects of the coupling between in-plane order and 
curvature (see, e.g., | p^ , p^ ). Here, in addition, the cou- 
pling between the nematic directors, m and n, and the 
curvature tensor K, of the form K : (m (g) m — n ® n) , 
produces interesting but complicated effects, which are 
out of the scope of this paper. In particular, shape fluctu- 
ations introduce an effective long-range coupling between 
director gradients [Q. On the other hand, the nematic 
tends to bend the membrane along its principal axes . 
Therefore the texture around a nematic bi-defect will de- 
form the membrane, and the membrane shape will react 
on the texture in a nontrivial way. 

The wall thickness can be estimated by assuming that 
the A term in Eq. (|l|) arises from anisotropic van der 
Waals interactions: A ~ A^.i'^ / {2TTd'^), where £ is the lin- 
ear size of the headgroup, d the membrane thickness, and 
Aa is the anisotropic Hamaker constant. Since Hamaker 
constants for interactions across a hydrocarbonic medium 
are of order T ||l5|, we take A^, ~ 0.1 T. Hence, with 
d ~ 40 A and £ ~ 10 A we find A ~ 2 lO"'^ Jm^^. Taking, 
e.g., ~ 3T, we obtain ^ = K'^/^ /{2\)^/^ ~ 1500 A. 
The wall thickness, which is of the order of 5i^, should be 
in the fim range. 

One way to produce flat nematic bilayers would be 
either to deposit the membrane on a water-air interface, 
or to compress a Langmuir monolayer of gemini until a 
second layer overlaps the first. Due to the micrometric 
thickness of the walls, striking defect patterns should be 
directly observable by optical microscopy. 

We thank A. Ajdari, P. Olmsted and D. Wu for useful 
discussions. LP acknowledges the support of a Chaire 
Joliot dc TESPCI. 



[4] F. N. Braun, T. J. Sluckin and E. Velasco, J. Phys.: 

Cond. Mat., 8, 2741 (1996). 
[5] M. R. Fisch and C. Rosenblatt, J. Phys. II France 4, 103 

(1994). 

[6] R. Oda, I. Hue and S. J. Candau, Chem. Commun. 2105 
(1997). 

[7] J.-B. Fournier, Phys. Rev. Lett. 76, 4436 (1996); J.-B. 
Fournier and P. Galatola, J. Phys. II France 7, 1509 
(1997). 

[8] S. J. Candau and F. C. MacKintosh, private communi- 
cation. 

[9] W. Helfrich, Phys. Rev. Lett. 21, 1518 (1968). 

[10] A similar decomposition has been introduced for a ne- 
matic subject to an inhomogeneous field by E. L. Aero, 
Crystallography Reports 40 823 (1995). 

[11] There are many equivalent field-line representations of 
the same texture, which are connected by uniform rota- 
tions of the directors. 

[12] F. C. MacKintosh and T. C. Lubensky, Phys. Rev. Lett. 
67, 1169 (1991). 

[13] T. C. Lubensky and J. Prost, J. Phys. II France 2, 371 
(1992). 

[14] We thank P. Olmsted for bringing this problem to our 
attention. 

[15] J. N. Israelachvili, Intermolecular and Surface Forces 
(Academic Press, London, 1992). 



Dipartimento di Scienze Fisiche and Unita INFM, Uni- 
versita "Federico 11", Mostra d'Oltremare, Pad. 19, I- 
80125 NapoU, Italy. Associato INFN, Sezione di Napoh. 

[1] P. G. de Gennes and J. Prost, The Physics of Liquid 
Crystals (Academic, New- York, 1993). 

[2] P. M. Chaikin and T. C. Lubensky, Principles of Con- 
densed Matter Physics (Cambridge U. P., Cambridge 
(UK), 1995). Chap. 9. 

[3] There is a wealth of experiments in thin nematic cells 
motivated by the development of nematic displays. See, 
e.g., S. V. Yablonskii et al., JETP Lett. 67, 409 (1998). 



4 



